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Abstract 

The Gilbert- Varshamov bound states that the maximum size A2 {n, d) 
of a binary code of length n and minimum distance d satisfies A2{n, d) > 
2"/V(7i, d — 1) where V{n, d) — X^iLo il) stands for the volume of a Ham- 
ming ball of radius d. Recently Jiang and Vardy showed that for binary 
non-linear codes this bound can be improved to 

for c a constant and d/n < 0.499. In this paper we show that certain 
asymptotic families of linear binary [n, n/2] random double circulant codes 
satisfy the same improved Gilbert- Varshamov bound. These results were 
partially presented at ISIT 2006 [3;. 

Index terms: Double circulant codes, Gilbert- Varshamov bound, linear 
codes, random coding. 



1 Introduction 

The Gilbert- Varshamov bound asserts that the maximum size Aq{n,d) of a 
q-aiy code of length n and minimum Hamming distance d satisfies 

AJn,d) > ■ . (1) 

' Eto(-)(9-i)^ 

This result is certainly one of the most well-known in coding theory, it was 
originally stated in 1952 by Gilbert [5] and improved by Varshamov in [15] . In 
1982 Tsfasman, Vladuts and Zink [14j improved the GV bound on the number 
of codewords by an exponential factor in the block length, but this spectacular 
result only holds for some classes of non-binary codes. Recently Jiang and 
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Vardy [6] improved the GV bound for non-linear binary codes by a linear factor 
in the block length n to 

2"' 

for d/n < 0.499, for a constant c that depends only on the ratio d/n and 
where V{n, d) = X]f=o (^) stands for the volume of a Hamming ball of radius d. 
This new bound asymptotically surpasses previous improvements of the binary 
Gilbert- Varshamov bound which only managed to multiply the right hand side 
in ([l]) by a constant (see [6j for references). The method used by Jiang and 
Vardy relies on a graph-theoretic framework and more specifically on locally 
sparse graphs which are used to yield families of non-linear codes (their result 
was later slighlty improved in [H]). In this paper we also improve on the the 
Gilbert- Varshamov bound by a linear factor in the block length but for linear 
codes, thereby solving one of the open problems of [6]. The method we use is 
not related to graph theory and relies on double circulant random codes. 

Double circulant codes are [2n, n] codes which are stable under the action 
of permutations composed of two circular permutations of order n acting si- 
multaneously on two differents halves of the coordinate set. These codes can 
also be seen as quasi-cyclic codes, a natural generalization of cyclic codes jl3] . 
Their study started in 1969 in [8] and since they gave some very good codes 
it was natural to wonder whether they could be made to satisfy the Gilbert- 
Varshamov bound. A first step in that direction was made by Chen, Peterson 
and Weldon in [IJ who prove that when 2 is a primitive root of the ring Z /pZ for 
p a prime, double circulant [2p,p] random codes satisfy the Gilbert- Varshamov 
bound; unfortunately it is still unknown (this is Artin's celebrated conjecture, 
1927) whether an infinity of such p exists. Later Kasami [9], building on this 
idea, extended the result of to the case of powers of such p, and obtained 
a bound which is worse than the Gilbert- Varshamov bound by an exponential 
factor in the block length (though a very small one). Later Kasami's work 
was generalized to other cases in [71 [TTt 112] . and, in particular in [2j, bounds 
were proven for certain classes of quasi-cyclic codes that are worse than the 
Gilbert- Varshamov bound only by a subexponential factor in the block length. 
In this paper, building anew on Kasami's idea we prove, by using a proba- 
bilistic approach, that randomly chosen double circulant codes not only satisfy 
the Gilbert- Varshamov bound with high probability, but also the same linear 
improvement as that of Jiang and Vardy ([2]). 

The paper is organized as follows: in Section [2] we cover the main ideas 
involved. We start by recalling the probabilistic method for deriving lower 
bounds on the minimum distance of linear codes (section l2.0p . then we introduce 
double circulant codes in section 12.11 and derive ([5]) an upperbound on the 
probability that a random double circulant code contains a non-zero vector of 
weight not more than a given w. In section [2^2] we study the probability that 
a given vector belongs to a randomly chosen double circulant code. Finally in 
section 12.31 we derive our improved lower bound on the minimum distance in 
the simple case when the codelength is 2p and 2 is a primitive root of TLjpTL: 
the result is given in Theorem [H 
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in Section [3l we develop our method in the more comphcated case of block- 
lengths 2p"^, p a "Kasami" prime, in order to obtain an infinite family of double 
circulant codes with an improved minimum distance. Section [3.11 starts by giv- 
ing an informal sketch of the content of section [3l which is intended to give some 
guidance to the reader and discuss the technical issues involved. Section 13.21 
shows how to derive our main result, which is Theorem [HI from a proposition 
on the weight distribution of a certain class of cyclic codes. Finally section [33] 
is devoted to a proof of this last proposition. 

Section [4] concludes by some comments and side results. 

2 Overview of the method, the simple cases 

2.0 The Gilbert Varshamov bound for linear codes and its im- 
provement 

To put the rest of the paper into perspective and introduce notation, let us 
recall how the probabilistic method derives the Gilbert Varshamov bound for 
linear codes. Rather than bounding the code size from below by a function of 
the minimum distance, as in ([2]), we fix a lower bound on the code rate and 
find a lower bound on the minimum distance. We limit ourselves to the rate 
1/2 case because it will be our main object of study. 

Let Crand be the random code of length 2n and dimension k > n obtained by 
choosing randomly and uniformly a n x 2n parity-check matrix in {0,1}"^^". 
The probability that a given nonzero vector x = (xi . . . X2n) is a codeword 
is clearly 1/2". Let w he a positive number, not necessarily an integer. We 
are interested in the random variable X{w) equal to the number of nonzero 
codewords of C^and of weight not more than w. In other words 

X{w)= Yl (3) 

xe-B2n('!«) 

where B2n{w) denotes the set of nonzero vectors x of V2n = {0, 1}^" of weight 
at most w, and Xy^ is the Bernoulli random variable equal to 1 if x G C^and 
and equal to zero otherwise. Now whenever we prove that the probability 
P(X(u)) > 0) is less than 1, we prove the existence of a [2n, fe, d] code with 
k > n and d > w. Since the variable X(w) is integer valued we have 

P{X{W)>0)<E[X{W)] = Yl E[Xx] = |i?2„H|P(xG Grand) 

= \B2n{w)\l^. 

Hence, for every positive integers n and w satisfying \B2ni'w)\ < 2" there exists 
a linear code of parameters [2n,n,d > w]. Reworded, we have the following 
lower bound on d, essentially equivalent to ([T|). 

Theorem 1 (GV bound) For every positive integer n there exists a linear 
code of parameters [2n,n,d] satisfying 

\B2nid)\ > 2^ 
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In the present paper we shall prove 



Theorem 2 There exists a positive constant b and an infinite sequence of in- 
tegers n and [2n, n, d] linear codes satisfying 

\B2nid)\ > 6n2". 

This result, equivalent to ([2]) for rate 1/2, will be obtained by again choosing 
random matrices, but from a restricted class, namely the set of parity-check 
matrices of double circulant codes. 

2.1 Double circulant codes 

A binary double circulant code is a [2n, n] linear code C with a parity-check 
matrix of the form H = [I„ | A] where is the n x n identity matrix and 
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a„_i . 
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There is a natural action of the group 'LjnTL on the space V2n = {0,1}^" of 
vectors x = {x\ . . . Xn, Xn+i ■ ■ ■ X2n) namely, 

Z/nZx^an ^ ^2n 
(j, x) ^ j • X 

where 

1 • X — ,X\... Xfi—\^ X2ni Xn+1 1 • • • 3^2n— l) 

and j • X = (j — 1) • (1 • x). The double circulant code C is clearly invariant 
under this group action. Consider now C to be the random code Crand obtained 
by choosing the vector a = (ao...an-i) randomly and uniformly in {0,1}"". 
As before, we are interested in the random variable X[w) defined by ([3]) and 
equal to the number of nonzero codewords of Crand of weight not more than 
w. We are interested in the maximum value of w for which we can claim 
that V{X{w) > 0) < 1, for this will prove the existence of codes of parameters 
[2n, n,d> w\. The core remark is now that, if y = j • x, then 

Xy = Xx 

where Xx (Xx) is the Bernoulli random variable equal to 1 if x G Crand (y ^ 
Crand) and equal to zero otherwise. Let now i?2n(^) be a set of representatives 
of the orbits of the elements of B2n{w)-, i-e. for any x S B2n{w)-, \{j • x,j G 
Z/nZ} n B'2n{w)\ = 1. We clearly have > if and only if X'{w) > 

where 

X'{w)= 
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Denote by £(x) the length (size) of the orbit of x, i.e. ^(x) = € Z/nZ}. 

We have 

By writing P{X{w) > 0) = F{X'{w) > 0) < E[X'{w)], together with (jlj) we 
obtain 

p(xH>o)<5: ^ (5) 

d|n wt(x)<to 
e(x)=d 

Suppose in particular that n is a prime, in that case orbits are of size 1 or n, 
and a w < n then clearly the orbit of x has size n for any x G B2n{w), so that 
([5]) becomes 

V{X{w) > 0) < E [X{w)] In. 
If we can manage to prove that 

E\X(w)\<\B2niw)\^ (6) 

for constant c, then we will have proved the existence of double circulant codes 
of parameters [2n, n,(i > w], for any w such that \B2n{w)\ < \n2'^. 

2.2 The behaviour of P(x G C^and) 

To prove equality ([6]) we need to study carefully the quantities E [X^ , for 
X G B2n{w), since 

E[X{w)]= Yl 

xeB2„{w) 

For X e V2n, let us write x = (x/^,x/j) with x/,,x/j £ {0,1}". Consider the 
syndrome function a 

f : V2n Vn 

X h->- o"(x) = x*H = crL(x) + (Tij(x) 

where o"2,(x) = x^, and cr/j(x) = x/j*A. 

For any binary vector of length n, u = (mq, • • • > '^n-i), denote by u(Z) = mo+ 
uiZ + ■ ■ ■ + Un-iZ"'~^ its polynomial representation in the ring F2[Z]/(Z" + 1). 
For any u E Vn, let C(u) denote the cyclic code of length n generated by the 
polynomial representation u(Z) of u. Since ctr^x) has polynomial representa- 
tion equal to x^(Z)a(Z), we obtain easily 

Lemma 3 The right syndrome ct^(x) of any given x G is uniformly dis- 
tributed in the cyclic code C{xji). Therefore, the probability P(x G Crand) that 
X is a codeword of the random code Crand 'is 

. P(X G Crand) = 1/|C(xh)| if XL G C{xr), 
• P(X G Crand) =0 if XL ^ C{xr). 
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2.3 The case n prime and 2 primitive modulo n 

If n is prime and 2 is primitive modulo n then, over F2[Z], the factorization of 
+ 1 into irreducible polynomials is 

+ 1 = (1 + Z){1 + Z + + . . . + z"-^) 

and there is only one non-trivial cyclic code of length n, namely the [n, n — 
1,2] even-weight code. Therefore V{X{w) > 0) = V{X'{w) > 0) < E[X'(u')] 
together with ([U and Lemma [3] give 

wt(x£)+wt(x/j)<'u; wt(Xj;,)+wt(x/j)<UI 
wt(xfl) odd wt(xfl) even 

wt(xi) even 

P(XH>0) < 2\B2n{w)\^. 

We therefore have the following result: 

Theorem A If p is prime and 2 is primitive modulo p, then there exist double 
circulant codes of parameters [2p, p,d > w] for any positive number w such that 

2\B2p{w)\ <p2P. 

Unfortunately, it is not known (though it is conjectured) whether there 
exists an infinite family of primes p for which 2 is primitive modulo p. Therefore, 
to obtain Theorem [5] we will envisage cases when n is non-prime. This will 
involve two technical difficulties, namely dealing with non-trivial divisors d of 
n in ([S]), and non-trivial cyclic codes C{xji) of length n in Lemma El 



3 An infinite family of double circulant codes 
3.1 Preview 

In this section we will study the behaviour of the minimum distance of random 
double circulant codes for the infinite sequences of blocklengths 2n introduced 
by Kasami : we will have n = p^ for suitably chosen p. We will first specialise 
inequality (0) to this case, for which all the possible orbit sizes i are powers 
of p, p^, s < m. Applying Lemma [3] will lead us to an upper bound ([13]) on 
F{X(w) > 0) that involves the weight distributions of the cyclic codes of length 
n. This upper bound can be essentially thought of as the same as ([7]), plus a 
number of parasite terms involving all vectors x = (x/,, x^) of B2n{w) for which 
both Xi and x/j are codewords of some cyclic code of length n that is neither 
the whole space {0, 1}" nor the [n, n — 1,2] even-weight subcode. The problem 
at hand is to control the parasite terms so that they do not pollute too much 
the main term i.e. the right hand side of ([7]). To do this, the crucial part will 
be to bound from above with enough precision terms of the form 

MC)A,{C)^ (8) 

i+j<w ' ' 
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where C is a cyclic code of length n and Ai{C) is the number of codewords of 
weight i. In section 13.21 we shall state such an upper bound, namely Proposi- 
tion O and show how it leads to the desired result which will be embodied by 
Theorem [8l 

Section 13.31 will then be devoted to proving Proposition O It is not easy 
in general to estimate the weight distribution of cyclic codes that don't have 
extra properties, but it turns out that for these particular code lengths of the 
form n = p"^, all cyclic codes C have a special degenerate structure. Either C 
consists of a collection of vectors of the form {x,x, . . . ,x) where x is a subvector 
of length n/p and is repeated p times, or C is the dual of such a code. Section [3^2] 
will have reduced the problem to the latter class of cyclic codes only. Ideally, 
we would like to claim that the cyclic codes C have a binomial distribution 
of weights, i.e. Ai{C) « 2^(^)) however this is not true, the cyclic codes C 
have many more low-weight codewords than would be dictated by the binomial 
distribution. The problem of the unbalanced couples (i, j), {i small and j large 
or vice versa) in the sum ([5]) is therefore dealt with by the trivial upper bound 
A.i{C) < (^) : Lemma [11] will show that these terms account for a sufficiently 
small fraction of \B2niw)\/2^ . Lemma [101 is the central result of section [3.31 
which gives a more refined upper bound on Ai(C) for i well enough separated 
from 0, i.e. i > ktl for constant positive k. Fortunately, we do not need Ai{C) 
to be too close to the binomial distribution, and the cruder upper bound of 
Lemma [TT] will suffice to derive Proposition [5l 



3.2 Reducing the problem to the study of the weight distribu- 
tion of certain cychc codes 

Following Kasami [9], let us consider n of the form n = p"^ where 2 is primi- 
tive modulo p and 2^"^ ^ 1 modp^. It will be implicit that all the primes p 
considered in the remainder of section [3] will satisfy this property. Let us also 
suppose m > 2, since the case m = 1 is covered by Theorem [H 

It is known [9] that the irreducible factors of -|- 1 in F2[Z] are 1 + Z 
together with all the polynomials of the form 

l + Q{Z) + Q{Zf + ---Q{Z)P-^ (9) 

for Q{Z) = Z,= ZP,ZP\..., ZP""'\ 

Since n is a prime power, ([5]) gets rewritten through Lemma [3] as: 

F(XM>0)<± E (10) 

s=l wt(x)<i« ^ ' ^ -"^l 
C(xi)cC(xfl) 

Note that x G V2n has orbit length £(x) < n if and only if both x^ and x/j 
are made up of p successive identical subvectors of length n/p. Equivalently x^ 
and xji each belong to the cyclic code generated by the polynomial 

Pn{Z) = 1 + + Z^^'/P + ■■■ + Z^p-^>Ip. (11) 
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Let C„ denote the set of those cycHc codes of length n whose generator poly- 
nomial is not a multiple of Pn{Z). All the other cyclic codes of length n are 
obtained by duplicating p times some cyclic code of length n/p. Therefore, for 
s = m, the inner sum in (jlOp can be bounded from above by: 

ceeni+j<w ' ' 



where Ai{C) denotes the number of codewords of C of weight i. Applying ([12 
recursively, we obtain from (jlOp 



771—1 ^ 

p{xH>o)<j2 E E (^') 



We now proceed to evaluate the righthandside of (jlSh . The most technical 
part of our proof of Theorem [2] is contained in the following Proposition. 

Proposition 5 There exist positive constants q, K, ci and 7 < 1 such that, for 
any n = p"^ with p > q, we have \B2ni'^Kn)\ < 2" and for any positive real 
number w, K < w/2n < 1/4, and for any cyclic code C of Qn, we have 

Ai{C)A^{C)— < ci — — 7 

i+j<w ' ' 

Suitable numerical values of the constants are q = 14^, K = 0.1, 7 = 1/2^/^, 
ci = 26/5. 

Before proving Proposition O let us derive the consequences on the proba- 
bility P{X(w) > 0). That will lead us to our main result, namely Theorem [HI 
the consequence of which is Theorem [3 We have: 

Lemma 6 There exists a constant C2 such that, for any n = p"^, p > q, and 
for any K < w/2n < 1/4, 

cee„j+i<«> ' ' 

A suitable numerical value for C2 is C2 = 4.3. 

Proof: From Proposition [5] it is enough to show that the sum "^q^q^ ^71-dimC 
is upperbounded by a constant for any 7 < 1. Choosing a code C in 
is equivalent to choosing its generator polynomial, and from the list ([9|) of 
irreducible factors of -|- 1, we see that if we order all possible generator 
polynomials by increasing degrees, we have 1 and 1 + Z, then 2 polynomials of 
degree at least p — 1, then 4 polynomials of degree at least p{p — 1), ... then 2* 



8 



polynomials of degree at least p[p — 1)* ^ and so on. Therefore, since n — dimC 
equals the degree of the generator polynomial, we obtain 

^ ^n-dimC < 1 + 7 + + ^ 2*7P(P"^)'~' 

CGCn i>2 



< 



i>2 

2 



< 1 + 7 + 27^-1+ (^] VjV 



< I + 7 + 27P-1 + 



p — 1 / ^-^^ 

^ i>i 

2 \ 7 



^p-i; (1-7)2' 

With the values 7 = 2^/5, ci = 2^/^ and p > 14^ given in Proposition we 
obtain that C2 = 4.3 is suitable. I 

From (jl3p and Lemma [6] we obtain that 

P(X(.) > 0) < + c. V (14) 

to deal with this last sum we invoke: 

Lemma 7 For any prime p > 14^ and for any positive number w such that 
\B2ri{w)\ < n2", we have 

p' \B2n/p-{w/p')\ ^ 2 

Proof: Choose p times a vector of length 2n/p and weight not more than w/p: 
concatenate the resulting vectors and one obtains a vector of length 2n and 
weight not more than w. Therefore |-B2n/p(^/p)P ^ \B2niu!)\ and we have 



P^ \B2n/Aw/pn\ ^ / |ij2n(u;)| \ ^/^^ ^ y.^ F ^ 

The result follows from routine computations. I 

We see therefore from (|14p and Lemma [7] that, if we choose w such that 
1-5271(^^)1 < 6n2"', for 6 < 1, then, provided the conditions of Proposition [5] are 
satisfied, we have P{X{vu) > 0) < bc2 + 2c2/p- For C2 = 4.3 and any p > 14^ this 
quantity is less than 1 when b < 0.23. The largest w for which \B2niw)\ < bn2"' 
is readily seen to satisfy K < < j which means that all conditions of 
Proposition [5] are satisfied, so that we have proved: 

Theorem 8 There exist positive constants b < 0.23 and q, such that for any 
prime p > q such that 2 is primitive modulo p and 1P~^ ^ 1 modp^, and 
for any power n = of p, there exist double circulant codes of parameters 
[2n,n,d > w] for any w such that \B2n{'w)\ < 6n2"'. A suitable value of q is 
q = 14^ and the first suitable prime p is p = 2789. 
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3.3 Proof of Proposition [5] 

Our remaining task is now to prove Proposition O We start by noting that 
Proposition [5] is stated with a positive real number w, because the discussion 
starting from (jl3p involves balls of non-integer radius. However, it clearly is 
enough to prove it only for integer values of w. 

The crucial part of the proof will be to bound from above the weight distri- 
bution of C, for C £ C„. Let us note that, since the polynomial Pn{Z) defined 
in (jlll) is an irreducible factor of + 1, the code C belongs to Cn if and only 
if Pn{Z) divides the generator polynomial of the dual code C"*-. This means 
that any codeword of must be obtained by repeating p times a subvector of 
length n/p. Equivalently, a generating matrix of C"*", i.e. a parity-check matrix 
of C is of the form 



meaning that it equals the concatenation of p identical copies of an r x n/p 
matrix A. 

We shall need the following lemma. 

Lemma 9 Let H^^ = [Ir | Ir I ■ ■ ■ I Ir] the r x tr matrix obtained by concate- 
nating t copies of the r x r identity matrix. Let atr be the associated syndrome 
function: 



Let w < tr be an integer. Then, for any s G {0, 1}^', the number of vectors of 
length tr and of weight w that map to s by atr is not more than: 



where w = totr. 

Proof: Let X be a random vector of length tr obtained by choosing indepen- 
dently each of its coordinates to equal 1 with probability uj. The probabilities 
that any given coordinate of atr(X.) equals or 1 are those of a sum of t inde- 
pendent Bernoulli random variables of parameter iv, namely: 



He = [A I A 



A] 



CTtr : {0, 1} 





1 + (1 - 2u;)* 
2 



and 



1 - (1 - 2ujY 
2 



Since all the coordinates of atr (X) are clearly independent 




(15) 



Now let W = wt (X) be the weight of X. We have 
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where h denotes the binary entropy function, h{x) = —x log2 x— {1 — x) log2(l — 
x). By a variant of Stirhng's formula [13j[Ch. 10, §11, Lemma 7] 



>2"^("V\/8na;(l-a;), (16) 



therefore: 

V{W = w)> ^ > ^ 



For given s, let A'^ denote the number of vectors of length tr and weight w that 
have syndrome s. Since F{atr(X.) = s \ W = w) = N^/ (JJ) we have 

P(at,(X) = s) > P(cJt,(X) = s\W = w) V{W = w) > 



) V2t^' 



Hence, by (fT5|) . 



, l + \l-2uj\^Y ftr 

< V2tr 

' ' w 



which is the claimed result. I 

Lemma 10 Let < k < 1/4. There exist q, such that for any p > q, n = p" 
and for any code C €z Qn, the following holds: 

• either C = {0, 1}" or C equals the even-weight code, 

• or, the weight distribution of C satisfies, for any i, ku <i < n/2, 

where r = n — dim C. 
For K = 0.07 a suitable value of q is q = 14^. 

Proof: If r = or r = 1, i.e. C equals the whole space {0, 1}" or the even- 
weight code, there is nothing to prove. Suppose therefore r > 1. From the 
factorization ([9]) of Z" + 1 into irreducible factors we see that we must have 
r > p — 1. From the discussion preceding Lemma [9] we must have 

r <n-p'^'^{p-l) =n/p (17) 

and a parity-check matrix of C is made up of p identical copies of some r x n/p 
matrix A. Therefore, after permuting coordinates, there exists a parity-check 
matrix of C of the form 

He = [B |Ir I Ir I • • • I Ir] 

where B is some r x (n—rt) matrix and is followed by t copies of the r x r identity 
matrix. The integer t can be chosen to take any value such that 1 < t < p: we 
shall impose the restriction 

t<p^/^. (18) 
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For any x G {0, 1}", write x = (xi,X2) where xi is the vector made up of the 
first n — tr coordinates of x and X2 consists of the remainding tr coordinates 
Now the syndrome function a associated to He takes the vector x G {0, 1}" 
to <t(x) = xi *S + (Ttr(x2) where atr is the function defined in Lemma [H The 
code C is the set of vectors x such that o"(x) = 0, therefore by partitioning the 
set of vectors of weight i into all possible values of xi we have, from Lemma (9) 



3 J J 



for any i such that 
Notice that: 

so that (fTOl) becomes 



i > tr. (20) 

tr \ (n-tr\ (!) {tr-j) ( r)> 



Ai{C) < V2t^ 




< V2tr{tr + 1) . max ' ' ' • (21) 



Set i = LTi and j = atr, we have: 



{tr)3'-{tr-j)\ 



< 



^ , . Lnn-1 since Q>in-trr- /{try. 



< 



{n-tryr^tr 



(1 _ «I 



tr\tr 



We have seen (I17p that r < n/p and t < p^^^ (condition (llSp ). therefore tr/n < 
l/p2/3 < 1/2: by using the inequality 1 — x > 2 valid whenever < x < 1/2, 
we therefore have 



(7) itr-j) 



and by using (*[) < 2*'^'*(") we finally get 



<22*'^'/^J(l-0*'--'(*j) 



0^ < 2t'-(^-s("lk)) 



KtrJ 
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where D{x\\y) = xlog2 | + (1 — x) log2 Together with ([2T]) we get: 



with 



/(i) = max 5(a, t) (22) 

0<a<l 

where ^(a, i) = log2(l + |1 - 2a|*) - t£>(a||i) (23) 

and P = I log2 \/2tr + ^ log2(tr + 1) + 2t^r/n. Write log2(ir + 1) < 1 + log2 tr 
to get /? < (| + I log2(ir))/r + 2t^r/n. By using t < p^/^ and p - 1 < r < n/p, 
we get 

3 log2 tr ^ 3 1og2(r + l)-^/3 ^ ^ log2(r + 1) ^ ^ log2P 
2r 2 r r ~ p — 1 

and 

3 21og2P 2 



2{p-l) p-1 pV3 ■ 

We see that /3 can be made arbitrarily small by increasing the value of p. A 
numerical computation gives us f3 < 0.152 for allp > 14^. 

Since we have supposed i < n/2, we have l < 1/2 so that the definition (j22p 
and (f23]) of / can be replaced by the equivalent 

/(O = max g{a,L) 

0<a<i, 

g{a, l) = log2(l + (1 - 2q)*) - tD{a\\L) 

from which we easily see that g and / are decreasing functions of l. We see that 
/(k) can be made arbitrarily small, for all k > 0, by choosing t big enough. 
Numerically, by choosing t = 14, k = 0.07 and p > 14^, we see that ()20p is 
satisfied and we get, for all 0.07 < i, f{i) < /(k) < 0.24. We obtain therefore 
that, for all nn < i < n/2, 

A,(C) <2-0-608-^'' 



which proves the lemma. I 

To prove Proposition [5l we need a final technical lemma, of a purely enu- 
merative nature. 



Lemma 11 Let < k < K < 1/4. There exist an integer uq and e > such 
that, for any n>nQ, w = 2ujn with K <uj < 1/4, 




For K = 0.07, K = 0.1, no = 14^, a suitable value of e is e = 0.004. 
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Proof: Clearly we have: 

^ J \J J \KnJ\w — Kn 



i+j<w 
i<Kn 



2n{2h{ui)-h{K)-h{2ui-K)) 



2n(2h{K)-h(K)~h{2K-K)) 



since 2h{oj) — — h{2u) — k) is an increasing function of uj. By (jl6p we have 
(w)|, SO that we obtain, since k < 1/4, 

' n\ f Tl\ ^ \B2n{w)\ ^\B2n{w) 

i+j<w 



2 V < ;.5/2 l"2nV«^;| 

I 7 / I - 2"(2/i(i<')-/i(K)-'i(2A'-K)) - 



i<Kn 



for any n > no with e < 2/i(i^) - - h{2K - k) - I 

Proof of Proposition [5) If C = {0, 1}" or if C is the even-weight subcode, 
then MO < (^), and T.^+,^^ MC)A,{C) < E.+,<. ( •) = I^SnMI- The 
result clearly holds for any ci > 2/7. 

Let C € Cn with r = n — dim C > 1. Let us write: 

^ A{C)A,{C) = S, + S2 

i+j<w 

with 

•^1=1^ E and ^2 = ^ ^ ^.(C)^,(C). 

Kn<i,j i<Kn 

By Lemma [TOl we have 

Q <^ V M M 1 ^ l-g2n(^;)| 1 

To upperbound ^2 we simply write Ai{C) < ("). By Lemma \TT\ we have 

^ |g2n(t^)| 2^ _ \B2n{w)\ 2' \B2n{w)\ 2^ 

2 - 2"- 2=" 2*^ (2£P)"/p ~ 2" (2^P)''' 

since we have seen pT|) that r < n/p. By choosing p > ^ we obtain 



S2 < 



5e 

\B2niw)\ 1 



2"- 2'"/5 ' 

This proves the result with 7 = 1/2^/^ and ci = 2^/^ 
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4 Comments 



The probabilistic method we used easily shows that almost all double circulant 
codes of the asymptotic family presented here satisfy an improved bound of 
the form ([2]). Actually we suspect that this is also the case for most choices 
of n : this is suggested by computer experiments with randomly chosen double 
circulant codes of small blocklengths. 

We have tried to strike a balance between giving readable proofs and deriv- 
ing a non-astronomical lower bound on the prime p in Theorem [8l In principle, 
the numerical values could be refined. In particular, the constant b of Theorem 
[8] could be made to approach 1/2 (as in Theorem but at the cost of a larger p. 
If we convert the formulation of Theorem [8] in the form ^ (which just involves 
switching from \B2nid)\ in Theorem [2] to \B2nid — 1)| in dJ])) we obtain a con- 
stant c which is of the same order of magnitude, but somewhat worse, than the 
improved constant c ~ 0.102 of pj6j for Jiang and Vardy's method. 

In this paper we only consider the binary case with codes of rate 1/2 but the 
method can be straightforwardly generalized to the case of different alphabets 
and to quasi-cyclic codes of any rational rate (though at the cost of a wors- 
ening of the constant b) by considering for parity check matrices vertical and 
horizontal concatenations of random circulant matrices. 

Finally, a natural question is to wonder whether the ideas developed in this 
paper can be extended to Euclidean lattices in a way similar to the generaliza- 
tion of Jiang and Vardy's method to sphere-packings of Euclidean spaces [TO] . 
A positive answer to this question is given in the paper [3]. 
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